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Abstract. We report the first results of the MUSIC project. It consists of two data 
sets of resimulates clusters extracted from two large dark matter only simulations: 
Marenostrum Universe and Multidark. In total, the MUSIC contains more than 400 
clusters resimulated with high resolution both with radiative and non-radiative physics 
included. Here we present the first results on the properties of the baryon content and 
the Sunyaev Zeldovich scaling relations. 



1. MUSIC Datasets 

The first dataset is composed of 164 resimulated clusters extracted from the MareNos- 
trum (MN) Universe. The MN Universe (Gottlober et al. 2006) is a non radiative SPH 
cosmological simulation of a 500 Mpc h~^^ cubic box with WMAPl cosmological pa- 
rameters {D.\= 0.7, Q.M= 0.3, Q.b = 0.045, h = 0.7, cr8=0.9 and a slope of « = 1 for 
the initial power spectrum). It contains a total of 2 x 1024^ particles equally divided 
between dark matter and gas. We selected for resi mulation a total of 82 relaxed clus- 
ters and 82 bullet-like clusters and using GADGET ("SpringellllOOl) ,with 8 times more 
particles {mDM - 1.03x10^ m^ and mg^s - 1.82xlO^/i~^ m©) and adding radia- 
tive physics (i.e. cooling, UV photoionization, star formation and galactic winds and 
SN thermal and kinetic feedbacks). The second dataset is composed of the 282 more 
massive (most of them with a mass M >l0^^h'^ m© at z = 0) galaxy clusters extracted 
from the MultiDark (MD) simulation. The MD simulation is an ART dark matter only 
simulation performed by Anatoly Klypin at NAS Ames, containing about 8.6 billion 
particles in a (1 Gpch~^)^ volume that was performed using WMAP5 cosmological pa- 
rameters (Qa= 0.73, Qm^ 0.27, Qt ^ 0.0469, h ^ 0.7, o-8=0.82 and n ^ 0.95). We 
resimulated them with 8 times more particles in a region of 6 Mpc centred around each 
object at z=0. We added SPH gas particles in the high resolution area and run them in- 
cluding radiative (as in the case of the MN dataset) and non-radiativ e physics. The mas s 
resolution is moM = 9.01x10^ h'^ m^ and m^^, = 1.8x10*^/j"^ ( Klypin et a"l1l200lh . 



In many of the resimulation areas we often found more than one cluster. Thus, the 
total number of resimulated clusters we finally compiled was over 500 with M > 10'"* 
m©. The two datasets have now joined to form the MUSIC (MareNostrum-MUltidark 
Simulation of galaxy Clusters) project. 
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2. Baryon properties in simulated galaxy clusters 

We analyzed some of the internal properties of the clusters. We mainly focused on 
the baryon properties, in order to see how relevant is the effect of overcooling on the 
clusters simulated with radiative physics. Overcooling is an hot topic of cluster simu- 
lations (Kr avtsov e t al. 2005): the huge star formation rate in the cluster center implies 
a star fraction larger than the observed one. We studied the behavior of a set of cluster 
properties (total mass, gas and star mass, baryon , star and gas fraction) at six different 
overdensities ranging from A = 200 to A = 2500, where A is the ratio between cluster 
density and local critical density. These two limit values correspond to an overdensity 
(A = 200) at which the cluster mass well approximates the virial mass and a value (A 
= 2500) often refeiTcd in X-ray observations allowing us to study the effect of over- 
cooling. We also analyzed the dataset at five different redshifts (from z = to z = 0.43) 
in order to study the effect of cluster evolution. We define the overdensity radius as 
the radius of the sphere within which the average density of the cluster is A times the 
critical density: 

-;rp,,(z)Ar^ = MrorirA) (1) 

The critical density is defined as pc(z)=3Ho^E(z)^/87: G depending on the assumed cos- 
mology. We calculated the gas fraction as f gas = MgasMroT', similarly, the star and 
baryon fraction as/,,,,,- = MstarMroT andfhar =(Mstar+Mf.as)/MTOT- Figures 1 and 2 
show the baryon, gas and star fractions along the redshift for clusters of the MD dataset 
simulated with radiative physics at two different o verdensities. In fa. i A = 500 and 
the mean values agree very well with observationsCLaRo que et all 2006 '): the baryon 



fraction is close to the critical ratio H^/Q,,, as expected and the gas fraction is about 
twice bigger than the star fraction, lnfig.2 A = 2500 the effect of overcooling (when the 
clusters are less evolved). The radiative clusters of MN dataset show a similar behavior 
but, being the masses smaller than those of MD dataset, the effect of overcooling is 
slightly stronger and the star fraction higher. On the other hand, morphology seems 
not to affect the baryon properties. Non-radiative clusters are obviously not affected by 
overcooling and the mean gas fraction approaches the critical ratio as we move towards 
A = 200. 



3. SZ scaling relations 

The i ntegrated brightness of the t hermal component of the Sunyaev-Zel'dovich effect 
(SZ) dSunyaev & Zeldovichlll97dh . Y, is one of the strongest proxy to estimate the mass 



of galaxy clusters.Y is the integration over the solid angle subtended by the cluster D. 
of the Compton-y parameter defined as the integration of the electronic pressure along 
the line of sight. We studied the Y-M scaling relation of our simulated cluster sample 
to check it and compare it with other simulations, theory and observations. Under the 
assumption of an isothermal model the integrated Y is proportional to the integral of the 
electron density Ue'. 

where is the electronic temperature and Da the angular distance If we assume that 
cluster formation is driven by a spherically gravitational processes mainly due to the 
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Figure 1 . Baryon (squares), gas (diamonds) and star (triangles) fraction vs redshift 
at overdensity A - 500 of MD clusters resimulated with non-adiabatic physics 
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Figure 2. Baryon (squares), gas (diamonds) and star (triangles) fraction vs redshift 
at overdensity A = 2500 of MD clusters resimulated with non-adiabatic physics 



dark m atter component with baryons and dark matter in h ydrostatic equilibrium (iKaiser 
19861) we find a scaling relation based on self-similarity (IBonamente et al.ll2006l) : 



YDi^fgasMZ'Eiz) 



(3) 



where E(z) is the scale evolution. If we study this scaling relation in the log space, in 
the form logfFj = A\og(X)+B we expect a slope A = 5/3 = 1.66. To estimate Y for each 
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cluster we built SZ maps by Compton-j parameter : 

y - y nejTejWQ r - rcdm I, (4) 

I 

where hi is the smoothing length of the particle and W the particle kernel, |r- Vcdm I the 
position of SPH the particle respect to the centre of mass. All Y-M scaling relations, 
independently from physics, redshift or overdensity, showed a slope at least very close 
to the expected value of A = 1.66 {seefig.3) 
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Figure 3. Y-M scaling relation of MD radiative clusters at z = 0, A = 500 

Concluding, we analyzed two datasets of clusters for a total of more than 600 
objects simulated with high resolution and both radiative and non-radiative physics. 
The effect of overcooling becomes relevant only at high overdensities and the Y-M 
scaling relations show a very good agreement with a self-similar model. 
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